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ABSTRACT 


We describe some techniques for the construction of three-dimensional composite overlapping grids 
using the grid construction program CMPGRD. The overlapping approach can be used to generate 
grids for regions of complicated geometry. The grids can be constructed to be smooth and free from 
coordinate singularities. The ability to create smooth grids for complicated regions is an important 
first step towards the accurate numerical solution of partial differential equations. We describe how 
to create grids for surfaces defined by cross-sections such as an airplane wing. We also describe how 
the patched surfaces generated from a CAD package can be used within the CMPGRD program and 
how grids can be created in regions where surfaces intersect. 


INTRODUCTION 


We describe some techniques for the construction of three-dimensional composite overlapping 
grids using the grid construction program CMPGRD. A composite overlapping grid consists of a 
set of logically rectangular component grids which cover a region and overlap where they meet. 
Interpolation conditions connect solution values defined on the grid. In the paper Composite Meshes 
for the Solution of Partial Differential Equations [l] we described CMPGRD and discussed techniques 
for the solution of elliptic and hyperbolic partial differential equations (PDEs) on overlapping grids. 
Our examples were limited to two-dimensional problems. In this paper we emphasize the extensions 
made to the code for the generation of three-dimensional grids. The overlapping approach can be 
used to generate grids for regions of complicated geometry. The region can be divided into a number 
of sub-domains for which a component grid can be more easily created. The component grids can 
be constructed to be smooth and free from coordinate singularities. Once all component grids have 
been created the CMPGRD program will automatically determine the interpolation conditions' which 
connect the grids. CMPGRD has a very general algorithm for doing this which supports any number of 
component grids overlapping in any order. CMPGRD can generate the overlapping grids appropriate 
for higher order interpolation, higher order discretization, cell-centred or cell-vertex grids and the 
sequence of grids which can be used for the multigrid algorithm [2]. Some other references which 
describe overlapping grids and the solution of PDEs thereon include [3] [4] [5] [6] [7], 

The problem of grid construction can be made simpler by using the composite overlapping grid 
technique. However, the task of creating component grids, especially in three-dimensions, can still 
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be very difficult. In this paper we discuss some ways to create three-dimensional component grids 
suitable for use with CMPGRD. The emphasis here is on techniques which are closely linked to the 
overlapping approach and to methods which avoid coordinate singularities. 

We first describe how to create grids for surfaces defined by cross-sections, such as an ellipsoid or 
wing, where the cross-sections may converge to a point at one or both ends. The basic underlying 
problem here is the creation of a composite grid about a sphere which is free from singularities. A 
solution to this problem is to cover the sphere with more than one patch. We use an orthographic 
projection to create grids over the north and south poles of the sphere. This approach for a sphere 
can be generalized to more general surfaces which are defined by cross-sections: a separate patch is 
placed on the end of the surface where the cross-sections converge to a point. 

A complicated object such as an airplane consists of multiple surfaces which connect to each 
other, such as when a wing joins a fuselage. Using the cross-section technique component grids can 
be created for the wing and fuselage but there remains the problem of connecting the wing to the 
fuselage. The connection may either be smooth, in which case a fillet grid is appropriate, or the 
connection may be a corner in which case the curve of intersection between the two surfaces should 
be an edge of a grid. We consider the latter case here. One way to create a grid for the wing which 
matches to the fuselage is to project the end face of the wing grid onto the surface of the fuselage. 
We present an example of this technique. We have developed another, more general, approach for 
creating grids in the region where surfaces intersect. This method begins by computing the curve of 
intersection between the surfaces and then reparameterizing the surface using the intersecting curves 
to define the portion of the surface to use. We give several examples to illustrate this approach. 

For the description of complicated three-dimensional regions there are many advantages to using 
a computer aided design (CAD) package. We show how the patched surfaces generated from a CAD 
package can be used within the CMPGRD program. Although the patched surface may be smooth, 
the parameterization of the surface may not be smooth and it is thus necessary to reparameterize the 
surface. These reparameterized surfaces can then be used in CMPGRD to create overlapping grids. 
We show an example of a grid for a wing connected to an engine nacelle by a pylon. 


CREATING THREE-DIMENSIONAL OVERLAPPING GRIDS 


A composite overlapping grid consists of a set of component grids. Each component grid covers 
a portion of the computational region. Component grids overlap where they meet. Solution values 
are matched by interpolation at the overlapping boundary. In order to determine how to interpolate 
between grids CMPGRD requires knowledge of the component grid mapping not only at grid points 
but also at all intermediate positions. Thus the component grid must be defined as a continuous 
mapping. 

Component grid mapping: For the purposes of CMPGRD a component grid is defined as a 
smooth mapping from a unit cube, r, to the computational domain x € R 3 . CMPGRD needs to 
evaluate this mapping at any point r in the unit cube and requires both the image of the mapping, 
x(r), and its derivatives <9x(r)/<9r: 

<9x 

Component Grid Mapping : r — * (x, — ). 

OT 

A component grid has a number of characteristics associated with it. For example, each face of the 
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component grid has a boundary condition. The face may be part of a true boundary, or the grid may 
be periodic, or the derivative of the grid may be periodic or the face may be used for interpolation. 

CMPGRD provides various features for creating two-dimensional and three-dimensional compo- 
nent grids. One approach is for the user to provide a complete description of the component grid 
which is supplied to CMPGRD as an externally defined grid (i.e. Fortran subroutine). In this case 
one must supply the mapping from the unit cube to the computational domain and the derivatives 
of this mapping. This approach is useful if a grid has been created by some other package. If the 
component grid mapping is only known at a set of grid points then the mapping can be defined 
everywhere by using interpolation. CMPGRD has such a component grid interpolation routine. 

Another way to define a three-dimensional grid is to first define a surface in three-dimensions as 
an externally defined curve. This surface can then be automatically extended in the normal direction 
to create a three-dimensional grid. Alternatively two offset surfaces can be blended to form a space 
filling grid. 

Sphere in a Box: As a first example of creating a three-dimensional composite grid we consider 
the problem of generating a smooth grid for the region exterior to a sphere and interior to a box in 
which the sphere sits. We describe two ways to create a composite grid for this problem. 

A grid around the sphere can be created using a mapping defined by the standard spherical-polar 
coordinate transformation: 


Sphere : (r ]y r 2) r 3 ) — ► {x u x 2y x 3 ) 

(0, <f>, R) = (27T7'j , n(r 2 - .5), R n + (R b - /i a )r 3 ) 

{x [ , x 2} £ 3 ) = (/fcos(0)cos(</>), R sin(0) cos(^), /?sin(^)) 


This is not the recommended method to create a grid since the mapping has singularities at the north 
and south poles. These singularities will spell trouble when one wants to solve a PDE on the grid. 
In fact CMPGRD must know which grids are singular in order for its algorithm to work properly. 

By using multiple patches it is possible to cover the surface of a sphere with grid transformations 
so that there are no coordinate singularities. At each pole we create a grid defined by an orthographic 
projection, see figure (1). Consider the set of points (s i, s 2} — R) on the plane which is tangent to the 
sphere (radius R) at the south pole. The orthographic projection associates a point (s 1 } s 2 , — R) on 
this plane with a point x = (.Tj , x 2 , £ 3 ) on the sphere as the intersection between the sphere and the 
line through the north pole arid the point (s 1 ) .S 2 , — /?). This point of intersection is given by 


where 


Orthographic : (sj, s 2 ) — > 

(*i,*a,*3) = 


{x u x 2 , -Tri) 

(pcos(0), psm(ft), -() 

, fv 2 si nf 2 -s 2 (o ' 2 — s 2 ) 

(o 2 + s 2 ) ’ (cv 2 + .<r 2 ) ’ (o ' 2 + .s 2 ) 


o = 2 R 


a 2 s 


P = a ^rr ~ 7 

cv- + s z 


l I l 

s, 4 s, 


c = a 


rv 2 — s 2 


a 2 4 s 2 


cos($) = — , sin (9) = — 
.s s 

C 2 4 p 2 = R 2 


To complete the definition of the grid transformation \vc define the relationship between the unit 
square coordinates and the ( s j , 52 ) coordinates, 

•S| = (r, - .5) , ,s 2 = (r 2 - .5) s h 
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where the constants s a , Sb determine the extent of the orthographic patch. For later reference we note 
the derivatives of the orthographic transformation: 


dp _ 

. P c 

1 ^ 
1 

P 2 

d C_ 

~P 


ds 

’ s n 

’ ds 

R 

’ dp 

C 


ds 

_ Si 

ds 

f2 

de _ 

sin(0) 

d9 cos(0) 

ds i 

s 

’ di 2 = 

s 

ds x 

s 

ds 2 s 


(1) 

( 2 ) 


The full three-dimensional component grid is defined by extending each surface in the normal 
direction. The box is covered with a simple rectangular grid. The composite grid is shown in 
figures (3) and (4). 

Surfaces defined by cross-sections: The method described in the previous section can be 
generalized to create grids for surfaces defined by cross sections. We assume that the surface is 
defined as a function of a periodic variable 0, 0 < 6 < 2tt, and an axial variable C, -1 <<< +1: 


x = f(0,C) 


When £ is fixed f(0, £ = Co) will define a cross-sectional curve on the surface. The cross-sections 
should converge to a point at C = ±1. We will also see that the cross-sections must also tend to an 
elliptical shape as ( ±1 in order that there be no singularity in the derivative of the mapping at 

the poles. 

The cross-sectional surface will be covered by three patches: a central cylindrical patch and two 
end patches defined by orthographic projections. The central patch is defined by the mapping 


C‘2 * (Pi, r 2) “ _+ 

X — 

. dx dx 

where 6 = 2n T\ , 


{x u x 2} x 3 ) 

m o 


(2tt 



df 

d< 


) 


C = C (r 2 - .5) 


The constant Cc < 2 determines the portion of the sphere covered by the central patch. When Cc “ 2 
the central patch will extend all the way to each pole. The orthographic patch at the ( = +1 pole is 
defined by 


C'l : (ri,r 2 ) -> ( x 1 ,x 2 ,x 3 ) 

x = f(0,C) 

dx 


uk p [■ . cos((9) di ^ 

lb, = 


dx_ 

dr 2 


= — | COS 

s 


idf sin(0) di i 

m pm } ~~r lp 9 ( } J 


where p ) s and 6 are the variables associated with the orthographic patch: 


S\ = (»’i - 5) s a , .s 2 = (r 2 - .5) s b , 5 2 = s 2 + s] 

as 


cos($) = — , sin($) = — 

s s 


P = 


a 2 + s 2 


C = or 


Of 2 — s 2 


a 2 + s 


2 > 


p l + c 2 = R 2 


(3) 

(4) 

(5) 

( 6 ) 
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The Jacobian derivatives of the mapping, equations (5), (6) are determined in the standard way 
using the chain rule 


<9x 
dr i 


d0_d f d(_ df 

dr i d6 dr i d( 

— sin(0) df —p 2 

s dO + ~R 


(a ,df 

^s(0)— . 


These derivatives are written in the special form of equations (5),(6). We assume that we have a 
subroutine which defines the cross-sections f (0, £) and also certain derivatives of the cross-sections: 


Cross-section routine: (0, £) 


^ 1 df df. 

p de ,p d<) 


( 7 ) 


We choose to have the routine return (J§|,/»§[) as opposed to (§J, §*) in order to avoid inaccuracy 
in numerically computing the Jacobian. In fact, in order that the derivatives of the patches exist at 
the poles the cross-sectional function f must satisfy 


lim^o 

lim J _, 0 


m df 

cos( % 


sin(0) df 
s dO 



r . ( a\ df , cos ( 0 ) m 


dO 


} =C 2 


( 8 ) 

( 9 ) 


where the constants C\, C \ must be independent of 9. These conditions will be satisfied provided 
that the cross-sections tend to an elliptical shape at the poles: 


/,(0, C) ~ p(acos(0) + 6sin(0)) . 


Condition (8), for example, follows from writing equations (5) in the form 


dx 

d^x 


and noting that 



sin(fl) dfi 

s del s 2 ^sR 


l)(scos(0))(s^) 


Um4(-4 - l) (s cos(0)) = 0 . 
p— o s l ' s K ' 


Ellipsoid: A grid for an ellipsoid with principal axes of lengths (a, b, c ) can be created using the 
cross-section function 


Ellipsoid : (0, C) — 1 ► (xj , x 2 , ^ 3 ) 

(xi,X 2 ,x 3 ) = {fufijs) = {apcos(d),bpsin{0),cQ 


with Jacobian derivatives defined by 

J Bh/m \Bhid» yu/de 

pd/,/d( pdf,/a< pdh/d<; 


— asin(0) bcos(0) 0 
— a(cos(9) —b(sin(9) cp 


The resulting grid is shown in figure (5). 



Wing with Joukowsky Cross-Section: A slightly more complicated problem involves the 
generation of a grid around an wing which sits in a box. For simplicity we define the wing to have 
cross-sections which are Joukowsky airfoils. At the wing tip the cross-sections are forced to become 
elliptical since this is a necessary requirement for the grid to be smooth. The mappings are thus 
defined as 


Joukowsky wing : (rj, r 2 ) — ♦ (xj, x 2 , x 3 ) 

(x!,x 2 ,x 3 ) = g{p) f 2 ( 8 , 0>9{p)f 3 ( 9 , 0 ) 

The cross section functions f 2 and f 3 are defined as Joukowsky airfoils: 

/2 + 1/3 = tv+l/w 

w = ae ,e 4- ide' 6 / 4 (p) 

h = g(p) - g'(0p 

The function g(p) (recall p = \Jl — C 2 ) is chosen to be g(p) = tanh(/?p) and causes the cross-sections 
to converge to a point. The function / 4 satisfies / 4 (p = 0) = 0 causing the airfoil to become elliptical 
at the tip and f' t (p = 1) = 0 causing the cross-section at (0/0 = (0, 1) to lie in the in the x 2 — x 3 
plane. We make use of this latter result when we attach the wing to a fuselage in the next section. 
For the airfoil shown the free parameters defining the airfoil were taken as 

a = . 85 , c/ — .15 , 8 = 15° = 15(7r/180) r , £ = 2 

In order to have sufficient grid lines in locations of highest curvature we stretch the grid lines at the 
leading and trailing edges. 

Wing on a fuselage - first approach: In this section we consider the problem of attaching a 
wing to a fuselage. Here we describe a relatively simple procedure to accomplish the joining. Later 
we describe another solution to this problem. The procedure is illustrated in figure (7) and consists 
of the steps: 

1. Create component grids for the wing and for the fuselage. 

2. Deform the end of the wing grid to be the same shape as the piece of the fuselage it will attach 
to. 

3. Position the deformed wing onto the fuselage. 

For step 1 we use the wing grid as defined in the previous section. For simplicity we define the fuselage 
to be part of a cylinder: 

Cylinder : (n, r 2 ) -*• (ii,x 2 ) 

9 = 7r (r] — .5) , R = f?i + R 2 r 2 

(x u x 2 ,x 3 ) = (/Zcos(0), y a -I- ybar 3 , Rsin( 6 )) . 

To perform step 2 we deform the end of the wing to lie on a cylinder. For our geometry this is 
accomplished with the mapping 

*1 <- *1 + - {Ri - \fR\ - A 
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which transforms the end of the airfoil onto the cylinder. The mapping has an exponentially small 
effect on parts of the wing which are far from the end (figure (7b). Finally we move the wing to join 
the fuselage. The resulting composite grid will have the topology shown in figure (7c). 

In general, to perform step 2 of deforming the end of the wing we will need to know how to map 
a plane (or the end face of the wing grid) onto the surface of the fuselage. It is also important to 
perform this deformation without changing the surface shape of the wing. The composite grid for 
a Joukowsky wing attached to a cylinder in a box is shown in figure (6). The figure shows a cut 
through the wing and cylinder. The interpolation points are shown. 

This example illustrates some important advantages and disadvantages of composite grids. The 
advantage was gained because the grid generation problem was simplified by decomposing the domain 
into simpler problems. The disadvantage came when the grids had to be joined together along a 
physical boundary. 

M6 wing: In this section we use the ideas of the previous section to create a grid for the ONERA 
M6 wing which is free from artificial coordinate singularities. The M6 wing does, however, have a 
sharp trailing edge. The geometry of the M6 wing is defined by a sequence of cross-sections. Each 
cross-section is defined by a set of points. We fit a cubic spline to each cross section, parameterized 
by pseudo- arclength. Thus if the cross-section at £ — C, is defined by the points 

(A(».<i)t/j(*,C)) : 

then we define the pseudo- arclength by S\ = 0 and 

Sj = 8j- 1 + Cy/(yj - Vj- 1 ) 2 + (Zj - ?>-i) 2 7 = 2, 3, ... , rij , 

where C is chosen so that s U) = 1. We fit cubic splines 

y,(s) = spline({(sj, ?/,)}) 

Z,{s) = spline({(sj, z 2 )}) 

We use (linear) interpolation between cross-sections to define the surface of the wing everywhere: 

f 2 (0,O = (l-a)F,(s) + «Ws) 

— (1 — a)Z,(s) + aZ t+i (s) 
s = 0/(2tt) 

cv = (C - C.)/(G+i - C.) , C. < C < C .41 

The cross-sections at the wing tip should converge to a point. Moreover in order that the tip be 
smooth the cross-sections should converge to an ellipse. To ensure these conditions we slightly deform 
the tip of the wing. The definition of the cross-sections for the wing are thus 

M6 wing : (rj, r 2 ) — ► {x l ,x 2) x 3 ) 

(xi,x 2 ,x 3 ) = (L(, tanh(/?/})/ 2 ($, C)> tanh(/?p)/ 3 (0, 0) 

where L is the length of the wing and where /? is chosen sufficiently large so as to not significantly 
change the shape of the wing tip. The composite grid for an M6 attached to a cylinder is shown in 
figure (8). 
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Component grids created from CAD Package: A common representation of surfaces in 
computer aided design (CAD) packages is that of a Coons patch [8], also known as transfinite in- 
terpolation [9]. The Coons patch is defined as a mapping from the unit square into three-space. In 
this representation the unit square (rj, r 2 ) is divided into an array of sub-patches. There are n t x 
such subpatches. Each subpatch is defined as a polynomial, mapping the unit square (iq, u 2 ) (of the 
subpatch) into x € /f 1 . 

77 u — 1 71 v — 1 

Subpatch: x u (?/ 1 ,u 2 )= ]T ]T ^n u ™ u " ( 10 ) 

771 — 0 n=t 

where x^ m are constants. The subpatch boundaries, rq = r j , , i = l,...,n t — 1 and r 2 = r 2j , 
j = 1, — 1, are equally spaced in the (rq, r 2 ) plane, 




Ut + 1 


The full mapping is defined as 


, n. 


r '2,t 


U 2 + 1 




x(r,,r 2 ) = 


x u( 


7'l.i + l — r l,i 


r 3- r 3j ) 

»*2 j + l “ r 2,> 


for 7 j , < rj < rj 1 + i 


and r 2j < r 2 < r 2j+1 


The smoothness of the multi-patch surface x(rj,7* 2 ) is obtained by placing constraints on the sub- 
patch mappings at the boundaries of the sub-patches. Often patches are not parametrically smooth 
but only geometrically smooth: the derivatives with respect to r i or 7' 2 may not be continuous across 
different subpatches. In general we require parametric smoothness of surfaces if we wish to use the 
grid to solve a PDE problem. In this case we may want to reparameterize the patch to create a 
smooth parameterization. We are working on a variety of methods for this purpose. In figure (9) we 
show a composite grid for some surfaces created from the CATIA computer aided design package. 
This grid shows a wing connected to an engine nacelle by a pylon. 


CREATING COMPONENT GRIDS FROM INTERSECTING SURFACES 


Curves of Intersection: The first step in defining a grid in the region where two surfaces 
intersect is to define and parameterize the curve of intersection. Even assuming that the curve of 
intersection is well defined, consisting of one connected component, the parameterization of this curve 
is not well defined. This parameterization will be important when we wish to create a grid. 

Suppose we have two surfaces, S t : [0, l] 2 — > R 3 , i = 1,2 mapping the unit square into three- 
dimensional space, 

s, = {x = S,(r,) : r, = (r t , s,) e [0, l] 2 } 


and suppose these surfaces intersect along some curve (figure (10)). Denote this curve of intersection 
by C : [0, 1] — ► R 3 , 

C = {x = C(s) : s e (0,1]}. 


As the curve C(s) traces out a path on the surface S, it will also trace out a path on the unit square 
which is mapped to the surface S,. Denote this two-dimensional curve on the unit square of S. by 
R.: [0,1] - [0,1 f. 


R, = {r = R,(s) : 5 e [0, 1]}. 
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The curves R,(s) satisfy 


C{s) = S,(R,(s)) s 6 [0, 1], *=1,2 (11) 

We see that the curve of intersection can be represented in three ways, C(s), Si(Ri(s)) and S 2 (R 2 (s)). 
For consistency we want a numerical representation for these curves so that the defining conditions 
(11) are satisfied precisely (up to round off errors). This means that when we evaluate the curve we 
must solve the nonlinear equations (11), using Newton’s method for example. That is given a value 
of s we determine (C(s), Ri (s), R 2 (s)} satisfying (11). We will keep an accurate approximation to 
the curve to use as a starting guess for Newton. One reason we require such a precise definition for 
the curve of intersection is that the curve may be used to define a grid which has very small spacings 
between grid lines. In this case a less accurate representation of C or R^ may lead to erroneous results 
when interpolation points are computed for the overlapping grid. 

One way to choose the parameterization of C(s) is to use one of rj, si, r 2 , s 2 to parameterize the 
curve. For example if the curve Ri(s) is a single valued function of Si, as in figure (10), then we may 
use s = rj. The equations defining C(s) would be 


C(s) = S,(r t ) i = l,2 

5 1 =5 


( 12 ) 


For given s there are 7 equations for the 7 unknowns C(s), = Ri(s), * = 1,2. These equations 

can be solved by bisection in the variable rj. Another way to choose the parameterization of C(s) is 
to base the parameterization on another curve which is close to C. Thus suppose we have another 
approximate intersection curve (^(s). Given s we define C(s) to be the point on C which intersects 
the plane which is normal to the curve at the point C/i(s), see figure (11). 

C(s) = S,(r.) i = l,2 

^(c( s )-c,( s )) = o 

These seven equations define the seven unknowns C(s), r, = R,(s), i = 1,2. 

Component grids defined using the curves of intersection: We consider a few examples 
showing how a component grid can be constructed making use of the curves of intersection. In 
these examples the curve of intersection C(s), between the two surfaces will become the edge of a 
component grid. The faces of the component grid adjacent to this edge will lie on the respective 
surfaces. To parameterize these faces we will make use of the curves R t (s), corresponding to C(s), 
which lie in the unit square coordinates. By construction, the grid we create will have faces that lie 
precisely (up to round off errors) on the surfaces. 

The geometry of the first example is shown in figure (12) where two concentric cylinders intersect 
two planes. The goal is to create the grid x = G(r) which has faces coming from the cylinders and 
planes as shown in figure (13). Figure (12) shows the 4 unit squares and the curves of intersection. 
To define the faces of G we use the portions of the surfaces which are bounded by the curves of 
intersection. In the unit square coordinates the faces will correspond to the regions bounded by the 
images of the curves of intersection. The mapping for the face, Gi(ri, 7 * 2 , 0) corresponding to the end 
at r 3 = 0 of G can be defined as the composite transformation, 


G 1 (r 1 ,r 2 ,0) = S 1 (P,(r ll r 2 )), 
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where the mapping P] in the plane can be defined, for example, by transfinite interpolation 

r = Pi(r!, r 2 ) = (1 - r 2 )Ri(r 1 ) + r 2 R^n). 

The other three faces are defined in a similar manner (figure 13). The grid G is then defined as a 
transfinite interpolation between these four faces 

G(ri, r 2 , r 3 ) = (1 - r^G^O, r 2 , r 3 ) + riG(l, r 2 , r 3 ) + (1 - r 3 )G 1 (r 1 , r 2 , 0) 4- r 3 G(rj, r 2 , 1) 

-{(1 - r,) [(1 - r 3 )G,(0, r 2 , 0) + r 3 G(0, r 2 , 1)] 

+ r i [(1 — r 3 )Gi (1, r 2 , 0) + r 3 G(l, r 2 , 1)] } 

Figure (14) shows some surfaces and a grid that have been created using these techniques. Note 
that this parameterization of the grid G will only be smooth if the parameterizations of the curves 
of intersection, C,(s) are defined in a consistent way. A major catastrophe occurs, for example, if 
opposite faces have parameterizations defined in opposite directions so the grid folds over on itself. In 
a less severe case the grid lines may become highly skewed. Even if the parameterization is still single 
valued we wish to avoid such cases since finite difference approximations are usually less accurate for 
skewed grids. 

Grid for two intersecting spheres: In the second example we create a grid to be used in the 
region where two spheres (actually spherical shells) intersect. This collar grid will be used to join the 
spheres and looks something like a torus. Two of the faces of this collar grid lie on the the surfaces of 
the spheres. The other two faces lie on spherical shells which are offset from the spheres. The collar 
grid is shown in figure (15). The composite grid for the two spheres in shown in figure (16). Only 
the grids on the surfaces of the spheres are shown but the collar grid can be recognized. 

Wing on a fuselage - approach 2: We use the techniques for creating grids from intersecting 
surfaces to connect a wing to a fuselage. We first define surfaces for a cylinder and a wing with 
Joukowsky cross-sections. Another wing-like surface is defined by offsetting the wing in the normal 
direction. The two wing surfaces are intersected with the cylindrical surface and a wing grid is defined 
in the region bounded by the intersections. The composite grid is shown in figure (17). 


CONCLUSIONS 


We have given a brief description of some techniques for the construction of overlapping grids in 
three space dimensions. We have shown how to create grids for surfaces that are defined by cross- 
sections. We have described some ways to create grids in the regions where surfaces intersect. These 
grids are defined by smooth transformations and are free from artificial coordinate singularities. The 
ability to create smooth grids for complicated regions is an important first step towards the accurate 
numerical solution of PDEs on such regions. The results presented in this paper show some of the 
potential advantages of using overlapping grids to create grids for three dimensional geometries. 
However, there is still much work to be done to make the grid construction problem easier and more 
automatic. 
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Figure 3: Overlapping grid for a sphere 



Figure 4: Overlapping grid for a sphere showing interpolation points 





Figure 5: Overlapping grid for an ellipse showing interpolation points 



Figure G: Overlapping grid for Joukowsky wing on a cylinder 
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Figure 7: Attaching a wing to a fuselage - approach l 



Figure 8: Overlapping grid for an M6 wing on a cylinder 




Figure 9: Overlapping g: 



Figure 10: Curve of in 
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Figure 17: Overlapping grid created using intersecting surfaces 
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